Omnivory can both enhance and dampen perturbations in food webs 



Jaroslav Ispolatov & Michael Doebeli 
Department of Zoology and Department of Mathematics, 



University of British Columbia, Vancouver B.C., Canada V6T 1Z4 
Corresponding author: Michael Doebeli, Department of Zoology, 6270 University 
Boulevard, Vancouver, BC V6T 1Z4, Canada, email: doebeli® zoology. ubc.ca 

Abstract 

We investigate how perturbations propagate up and down a food chain with and without self-interaction 
and omnivory. A source of perturbation is a shift in death rate of a trophic level, and the measure of pertur- 
bation is the difference between the perturbed and unperturbed steady state populations. For Lotka-Volterra 
food chains with linear functional response, we show analytically that both intraspecific competition and 
intraguild predation can either dampen or enhance the propagation of perturbations, thus stabilizing or 
destabilizing the food web. The direction of the effect depend on the position of the source of perturbation, 
as well as on the position of the additional competitive and predatory links . These conclusions are con- 
firmed numerically for a food chain with more realistic Type-II functional response. Our results support the 
positions of both sides in the long-standing debate on the effect of intraspecific competition and omnivory 
on the stability of trophic systems. 
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I. INTRODUCTION 



flfl 



Food web stability is a topic with a tremendous scientific and practical significance 111, |3|, |4, 



13L 1141 11511 . A large number of both theoretical and experimental studies have addressed 



the question of how the stability of a food web depends on its structure, and specifically, on the 



presence of omnivory (see e 



.g E, HQ HQ, 



14l 1 171 . Il8h . just to mention a few). Omnivory, 



defined broadly as feeding on more than one trophic leve 



is a common phenomenon in natural communities lllq . 



, which generally includes cannibalism, 



1711 . However, a consensus on the effect of 



omnivory on food web stability is still lacking. Thus, according to some results from mathematical 



food web theory, omnivory destabilizes ecological communities [|9|,LL4|], whereas other models and 
experiments suggest that omnivory should be a strongly stabilizing factor in food webs I5I. Isi fioi 



120 • The notion of "stability" itself is broadly understood: the most common examples include a 



scarcity of secondary extinction events in response to removal or a shift in birth or death rate of 
a species [6] , temporal characteristics of a return to a new equilibrium following such removal 
[|5j], the mere existence of locally stable steady states 112]], and a global permanency of food web 



dynamics II 111 . 



In addition to this apparent lack of consensus on how food web structure affects stability, it 
seems that there still exist fundamental questions related to the effect of omnivory on food web 
stability which have remained unanswered, and even more, not specifically formulated. For exam- 
ple, in many practical situations a population at a certain position in a food web is harvested on a 
fairly permanent basis (for example, by fishing). Alternatively, reproductive rates or death rates of 
some species may shift due to environmental changes, endemics, etc. Such long-term removal (or 
possibly, addition, |20) of a resource from a certain trophic level naturally affects populations of 
the other levels in the food web. Generally, such removal may lead to extinction of certain species 
and to destabilization of initially steady state population dynamics in the food web. However, 
if the initial or unperturbed network is sufficiently stable, a continuous harvesting within certain 
limits may often result, after some transient behaviour, in constant shifts in population of other 
trophic levels, bringing the whole food web to a new steady state. The susceptibility of the steady 
state food web populations, or the difference between the new and old, unperturbed steady states, 
is another characteristic of the stability of food webs with respect to shifts in death or birth rates 
or harvesting. This new characteristic has an evident practical significance. For example, it is 
interesting and important to determine "sustainable" levels of harvesting that produce a limited 
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effect on the rest of the food web and on the harvested level itself. Intuitively it is clear that such a 
susceptibility-based definition of stability is related to previously used definitions, because a food 
web that is less susceptible to harvesting is less likely to leave a basin of attraction of a fixed point, 
a permanence domain, or suffer an extinction. 

As with other definitions of stability, susceptibility of a food web to long-term removal of 
resource from a certain trophic level depends on the position of this level in the web, predation 
strengths and natural death rates. The induced population shifts also depend on the structure 
(topology) of the food web, and specifically, on the presence, position, and strength of omnivory 
interactions. This leads to the formulation of the main question we address in this paper: How 
is the effect of the steady removal of a population of a certain level on the rest of the foodweb 
mediated by the presence of omnivory and self-interaction? In other words, we compare the 
steady state response of two model food webs, with and without omnivory trophic relations, to a 
shift of a death rate of a certain trophic level. When the response of the food web in the presence 
of omnivory is less pronounced than without it, omnivory is considered to be a stabilizing factor. 
Conversely, when the induced changes in population densities s are greater with omnivory than 
without, omnivory is considered destabilizing. To analyze this question in its most isolated and 
fundamental form, we focus on the topologically simplest form of a food web, a linear food chain. 
Also, our principle analysis deals with the most basic form of predator-prey interactions, i.e., 
with Lotka-Volterra models with linear functional responses, which turn out to be analytically 
tractable. Harvesting of a certain species is modeled as an introduction of a certain per capita rate 
of removal, which effectively modifies the death rate of the corresponding species in the food web.. 
We consider two possible forms of omnivory for linear food chains. First, we consider competition, 
or a negative self-interaction, at a certain trophic level, which is modeled as an additional per 
capita death rate term proportional to the population density at this level. This can be thought of as 
omnivory because it is mathematically equivalent to cannibalism with a linear functional response. 
Second, we consider intraguild predation, which we define as a "shortcut" trophic link introducing 
an additional predator-prey interaction between two non- adjacent trophic levels. 

Our conclusions confirm both existing viewpoints regarding stability and food web structure 
(and hence support neither as an overarching principle): both intraspecific competition and in- 
traguild predation can either stabilize or destabilize a food chain. Whether stabilization or desta- 
bilization occurs depends on the relative position of the level at which the food web is perturbed 
on the one hand, and where the additional interaction links are introduced on the other hand. We 
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also show numerically that both stabilizing and destabilizing effects of omnivory, derived analyti- 
cally for Lotka-Volterra systems, qualitatively hold for food chains with more realistic, non-linear 
functional responses. 

The paper is organized as follows: In the next section we define the framework of our food web 
study, which is based on systems of coupled Lotka-Volterra equations, and find the steady state 
solution for the population densities at all levels of the food chain. Next we consider the effect on 
the steady state induced by a shift of a death rate of a given trophic level and determine how this 
perturbation spreads over the food chain. We then derive how competition, or self-interaction, at 
a certain trophic level, modifies the steady state of a food chain, and consider the effect of such 
self-interaction on the spread of perturbation induced by a shift in a death rate. We show that 
depending on the position of self-interacting and death-rate modified levels, self-interaction may 
inhibit or enhance such propagation. Then the same analysis is repeated for omnivory links, and 
we present two examples of how such a link inhibits and enhances spread of perturbation. In the 
discussion, we present examples of quantitative similarities between the effects of omnivory and 
self-interactions in food chains with type-II functional responses food chain and our results for the 
Lotka-Volterra systems. We also illustrate how our results can be applied to more topologically 
complex food webs. The Appendix contains the description of a model with type-II functional 
responses in the presence of omnivory links. 



II. LINEAR LOTKA-VOLTERRA CHAIN 



We consider a linear food chain with nearest neighbour predator-prey interactions of Lotka- 
Volterra form. We do not limit ourselves to a specific number of levels in the chain and label the 
basal level by and the highest predator by n. The rate of change of the population density Xk(t) 
of trophic level k (which will also be called "species" k) is given by 

—j^- = Xk [^kdk,k-lXk-l — Gfc+l.fc^fc+l — dk] ■ (1) 

Here a i: j determines the strength of predation of species i on species j, \ k is the conversion 
efficiency which connects the birth rate of species k to the amount of other it consumes, and <4 is 
the per capita death rate of species k. Evidently, the top predator, occupying the nth level, does 
not have any species that prey on it, which could be expressed by setting x n+ i = 0. In addition, 
the basal species, occupying the lowest (zeroes) trophic level, is characterized by a logistic growth 
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term, which mimics the finite input of energy into the system, given in the form of a finite carrying 
capacity K for the basal species. Without loss of generality, the linear death term for the basal 
species can be absorbed into the linear part of the birth term j3, which yields the following for the 
rate of change of the population density of the basal species: 

dx 



dt X ° 



^ i 1 ~ ?) ~ 01,0X1 



(2) 



We look for a steady state solution of this system, dx* /dt = 0, which is defined by the following 
system of linear equations: 

* * Qfc+i,fc dk 

X k-1 — X k+1\ r - {J) 

AfcOfc^-i ^k a k,k-l 

Recurrence relations of the form x* k = F(x* k+2 ) given by Eqs. © indicate that in the steady 
state the Lotka-Volterra food chain splits into two subsets or partitions: the steady state of even- 
numbered levels depend only on another even-number level and a set of constants, and similarly, 
odd-number levels are coupled only to odd-number levels. To solve the system © one needs to 
consider the closing equations which describe the steady state populations of the top and bottom 
levels. Because the dynamics of the top predator population depends only on the first closing 
equation (or the boundary condition for the recurrence relation ©) defines the population of the 
n — 1 level, 

<-> = T~7T — ' (4) 

^n""n,,n— 1 

Consequently, sequential application of recurrence © allows one to determine the populations of 
lower levels of the same parity as n — 1 , down to level one or zero. The second boundary condition, 



«i,o V K 

which follows from the stationarity of the basal population, connects the populations of levels one 
and zero, thus linking the even and odd partitions. Finally, applying © upstream, or expressing 
x* k+1 through x^_ v we determine the steady state populations of levels of the same parity as n. 
Thus Eqs. © @J 0) define the stationary solution of the whole Lotka-Volterra food chain: The 
recurrence starts from the n — 1th level, descends to level one or zero depending on the parity of 
n, and then ascends to the top predator level n. This is illustrated schematically in Fig.Q] 

In the following we will widely use this natural separation of the food chain into two partitions: 
The regression in the "descending", or ll {n — 1} " partition starts from the n — 1st level and goes 
down along the indexes of the same parity as n — 1 to level one or zero, depending on whether n 



N^*=d 4 ^4 a 43 



• — e — • e- 



FIG. 1 : Direction of recurrence for finding the steady state populations x* k in a 5-level (0 - 4) Lotka-Volterra 
linear chain. The solution starts with finding the population of the second from the top level, x% and is 
propagated downwards so that x\ is determined via ((3]>. Then Xq is determined using the Eq. ©, and the 
recursion is propagated upward to determine xlj and x\. Here and below the levels of {n} ( or ascending) 
partition are shown as black circles and the levels belonging to the {n — 1} (or descending) partition are 
shown as white circles. 

is even or odd. It is important to note that the steady state population of a level of the descending 
partition is determined only by the level immediately above it in the same partition. Conversely, 
the regression in the"ascending", or ll {n} " partition, which spans indexes of the same parity as n, 
starts from zero or one and goes up to n. The population density of a given level in the ascending 
partition depends only on the population of its lower nearest neighbour in the same partition. 

Assuming for simplicity equality of the predation intensities, = a, conversion coefficients, 
Aj = A, and death rates, di = d, it is possible to express the general solution for the steady 
state populations in a compact form. The system of equations for the steady state populations © 
becomes 

A4_! = x* k+1 + -. (6) 
a 

The notation 

y k = x* k \- k / 2 a/d, (7) 

reduces the recurrence relation © to 



y k -i = y k+ i 



Thus, y k can be expressed as a sum of a geometric progression, and the general solution for 
arbitrary i and k reads 

4+2* = AX - - (8) 

CL 1 — A 

The first boundary condition defines the population of the n — 1 level, 

* d 

Xn ~ l = AV (9) 
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and, via ®, the populations of the levels of the {n — 1} partition. The second boundary condition 
links {n} to {n — 1} partitions, 



and yields the full explicit form of the stationary state populations. For example, for odd n, 



In the following we denote, where possible, the levels of the {n — 1} partition as % and levels of 
the {n} partition as j. It also follows from (fTT|) that for a given set of rate constants, the maximum 
length of the chain is finite and limited by the requirements that the descending partition does 
not exceed the carrying capacity even at its maximum, x* < K and the ascending partition does 
remain positive even at its minimum, x* n > 0. 

III. PERTURBATIONS CAUSED BY SHIFTS IN THE DEATH RATE OF A SINGLE LEVEL 

Now we consider how changes of conditions at a certain trophic level affect the steady state 
populations of all levels of the Lotka-Volterra food chain. We assume the simplest form of a pertur- 
bation equivalent to imposing a certain probability per unit time ea for each individual inhabiting 
the level p to be removed from the population. This can be taken into account by modifying the 
death rate coefficient, d, ' = d p + ea (the factor a is introduced for further convenience). 

We now investigate how such a shift in death rate of the level p affects the steady state popu- 
lation of all levels, or, in other words, how the perturbed steady state populations x' k are different 
from the unperturbed ones x* k . For simplicity in the following we will assume that predation in- 
tensities and conversion coefficients are constant for all levels, a^ = a and Aj = A. Evidently, the 
shift in a death rate affects only the levels in the direction of recurrence from the perturbed site p. 
That is, if a level of the {n} partition is perturbed (so that p E {n}), the perturbation first manifests 
itself at the p — 1 level, while x* p+1 and populations of higher levels of the {n — 1} partition remain 
unaffected. It follows from Eq. © that 




(10) 



d A fc -(™ +1 )/ 2 - 1 
a 1 — A 




(11) 



Xxl 



■p-i 



= x. 



* 




(12) 



and 




(13) 
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Consequently, the perturbation propagates down the food chain to the zeroth level if n is odd and 
to first level if n is even, 

A = x*, i>p; x\ = x* + A(p+ i_ t)/2 , i < P- (14) 
Then at the basal level, due to the second boundary condition ©, the perturbation changes sign, 

x ' x = x '~ ^RW^ (15) 

for odd p and n, and 

for even p and n. For the {n} partition the "ascending" regression yields 

x j = x j ~ 7^ A( p_ j+1)/2 ( 17 ) 
Here we introduced a universal notation for odd and even n, 

^db? ; <18) 

with ip — 1 for odd n and ip — — 1 for even n. Thus, when the death rate of a species in {n} 
partition is increased, the population of the levels of {n — 1} partition lying below p increases, 
while the population of all levels of the {n} partition, including the perturbed level p, decreases. 

Similarly, when the death rate of the level p E {n — 1} is shifted, the perturbation propagates 
recurrently upwards from j — p + 1, 

x '. = x l x >. = x*, j < p, x] = x* - e\ k , j = p + 1 + 2k. (19) 

In this case the {n — 1 } partition, including the directly affected level p remains unperturbed, while 
the population of the levels of {n}n partition that are above p is decreased. 

We illustrate the propagation of perturbations in a linear Lotka-Volterra food chain in Figs. [2J 

El 



IV. SELF-INTERACTION 

Self-interaction at level s in a linear food chain is taken into account by adding —aax 2 s to the 
right-hand side of the rate equation for x s , 

= x s [Aax s _i — ax s+ i — d — aax s ] . (20) 
8 



FIG. 2: Propagation of perturbation caused by a shift by ea of the death rate of a level belonging to the {n} 
partition; the recurrence relation (fT2l . perturbed by the shift in death rate, is shown by a dashed arc. The 
resulting spread of perturbation is shown by solid line arcs. The numbers above and below vertices indicate 
the change in the population of the corresponding trophic level induced by the perturbation. 

o 








FIG. 3: Propagation of perturbation caused by a shift by ea of the death rate of a level belonging to the 
{n — 1} partition; notations and symbols are the same as in Fig. [2j 

Here a defines the intensity of self-interaction, relating it to the intensity of predation a. As 
mentioned, the self-interaction can be thought of as either competition or cannibalism. As a con- 
sequence of self-interaction, the sth equation for the steady state population is modified to 



d 

Ax s _i = x s+ i H h ax s . 

a 



(21) 



Thus, similarly to the case of death rate shift, the self-interaction of a level s of one partition affects 
the population of those levels of the opposite partition which are recursively downstream from the 
perturbation site, that is, s — 1 if s E {n} and s + lifs E {n — 1}. We denote 5xk the change 
to the steady state population induced by a self-interaction, so that the new populations in system 
with self-interaction are x% + Sx^. Then, from (1201) , 



XSxs^i = 5x s+ i + a(x* s + 5x s ). 

The unperturbed steady state populations x* k are given by CCD- 
There are two principal cases: 

• s E {n}. Here self-interaction directly affects the s — 1 level, 

A5x s _i = + 8x s ), 



(22) 



(23) 



and recursively affects all lower levels i < s — 1 of the {n — 1} partition, 



\5xi-2 = Sxi (24) 

and all levels j of the {n} partition, 

\5%j = 5xj + 2 (25) 
From the boundary condition (flOl) and regressions (l24l |25T) it follows that 

5x s = - 7 ^v / A5s s _ 1 . (26) 

Thus 

tei = rr (2V) 

for i < s, and 

Q/X* 

Sx J = -l^T, u,'^ v, (28) 

for any j. We denote by a* = a^/VX a frequently occurring group of parameters which 
is proportional to a. 

• s E {n — 1}. Here self-interaction directly affects the s + 1 level and recursively affects 
only the higher levels of the same {n} partition. Hence 

5xj = -a\ {j ~ s - 1)/2 x* s , j > s (29) 

Thus a self-interaction in {n} partition increases the population of the lower levels of the {n — 
1} partition but decreases the population of all levels in {n} partition. A self-interaction in {n — 1} 
partition decreases the population of the higher levels of {n} partition and does not affect the 
{n — 1} partition. 



V. DEATH RATE SHIFTS WITH SELF-INTERACTIONS 



Having established how perturbations induced by a shift in a death rate propagate in the linear 
Lotka-Volterra chain and how the linear chain is affected by self-interaction, we now consider the 
first of our two principal questions: How does self-interaction affect the propagation of perturba- 
tions? Specifically, we would like to know if self-interaction stabilizes the system by reducing 
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FIG. 4: The effect of self-interaction at level 1 (shown by a large dashed circle) on propagation of perturba- 
tion caused by a shift by ea of the death rate at level 3; both self-interacting and perturbed levels belong to 
the {n} partition. The recurrence relation modified by self-interaction (1211) is shown by a dotted arc 

the amplitude of perturbations induced by a shift of death rate, or, on the contrary, whether self- 
interaction destabilizes the system by increasing this amplitude. We denote by Ax k the change in 
the steady state population of level k resulting from the death rate shift d' = d + ea in a chain 
with self-interaction of the form given by the Eq. (1201) at the level s . We are interested only in the 
terms that vanish when the death rate shift is zero, thus Ax k does not contain e-independent terms 
proportional to x*. Formally, the recurrent relation for Ax k follows from Eqs. (fl~2l 



AAx fc _i = Ax k+1 + 5 ktP e + 5 k;S aAx k , (30) 

where Kronecker's delta-symbol 5 iy j = 1 for i — j and 5 iy j = when i ^ j. Since death rate shifts 
and self-interaction can occur at levels belonging to any of the two partitions, there are 4 classes of 
scenarios. It turns out that among these 4 scenarios there are two examples of a stabilizing effect 
of a self-interaction, one of a destabilizing effect, and one of an absence of any effect on stability. 

• Both the perturbed level p and the self-interaction level s are located in the ascending {n} 
partition, Fig.Hl First consider the case when p > s. For i > s the levels experience only the 
effect of the death rate perturbation and Ax» is given by (fl~4~l) . For the s — 1st level Eq. (l30l) 
reads 

\Ax s -! = Ax s+1 + aAx s , (31) 

where Ax s+1 = e/X^ 8 ^ 2 from d). Similarly to (|27l[28|). we express Ax s through Ax s _i 
and obtain for i < s 



AXl ~ (l + a^A^ 1 -)/ 2 ' (32) 



and 



AXj ^ n + a*U(p+!-i)/ 2 (33) 



(1 + c**)A(p +1 -J)/ 2 

11 
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FIG. 5: The effect of self-interaction at level 2 of the {n — 1} partition with a perturbation caused by a shift 
of the death rate at level 3 of the {n} partition. 

for any j. Now we consider the case when p < s. Then self-interaction enters the recurrence 
relation first, 

AAx s _i = aAx s . (34) 



The equation for the level with a shifted death rate becomes 

aAx K 



\Ax 



p-i 



(35) 



\(s-p)/2 + e 

Expressing Ax s through Ai p _i using the boundary condition at the base of the chain, we 
obtain 

(l + a*)A 

Consequently, the results for i < p and all j are given by Eqs. (1321331) . Thus we ob- 
serve that when both the perturbed and self-interaction levels belong to the {n} partition, 
self-interaction weakens the effect of death rate shift on other levels, and therefore has a 
stabilizing effect. 

The perturbed level p belongs to the ascending {n} partition and the level of self-interaction 
s belongs to the descending {n — 1} partition, Fig. [5] If p < s, no joint effect occurs as the 
self-interacting level remains unperturbed by the death rate shift. However, in the opposite 

case of p > s, Ax { = e/A (p+1 ^ )/2 , i < p., see d), and Axj = -7^e/A (p - J+1)/2 , j < s, 
see (fT7l) . For j > s the effect of self-interaction simply adds to the effect of perturbation, 



Axi 



"7 



1 + 



a: 



^(p-j + l)/2 I rylf>\l/2 J ■ 



(37) 
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FIG. 6: The effect of self-interaction at level 3 of the {n} partition when perturbation is caused by a shift of 
the death rate at level 2 of the {n — 1} partition. 

Here we observe that self-interaction enhances the effect of perturbation, and hence is desta- 
bilizing. 

• The perturbed level p belongs to the {n — 1} partition and the level of self-interaction s 
belongs to the {n} partition, Fig. [6] There is no joint effect when p > s. When p < s, 
the correction to the j > p levels of {n} partition, including the level s, comes both self- 
interaction and shift of death rate. Hence, similarly to (1271) , we write for the i < s levels of 
{n — 1} partition 



Ax,; 



ae 



1 + a* 



and, similarly to (|28T) . for j < p levels of {n} partition 



ae 



+ a* 



_^Q-p)/2-l_ 



For j > p, a direct correction from the death rate shift is added, 



ae 



l + a* 



(38) 



(39) 



(40) 



Thus the self-interaction dampens the propagation of perturbations and stabilizes the system. 

• Finally, both the perturbed level p and self-interaction level s are located in the descending 
({n — 1} ) partition, Fig. |7J In this case no joint effect of perturbation and self-interaction 
occurs since the descending partition levels remain unperturbed and do not add any terms 
depending on e through self-interaction. 

A qualitative summary of the effect of self-interaction on propagation of perturbations is pre- 
sented in the Table 1 . 
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FIG. 7: The effect of self-interaction at level 3 of the {n} partition when perturbation is caused by a shift of 
the death rate at level 2 of the {n — 1} partition.. 





s G {n} 


s £ {n - 1} 


p e {n} 


dampens 
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p G {n — 1} 


dampens 


no effect 



VI. OMNIVORY SHORTCUTS 

Here we address the second principal question of our study: How does an intraguild trophic re- 
lation, or a predator-prey interaction connecting two non-adjacent trophic levels, affect the spread 
of perturbations in the food chain (where we again consider perturbation induced by shifting the 
p-level death rate)? A shortcut trophic link connects what we call an upper level u to a lower level 
w, u — w > 2. As in the self-interaction case, we would like to know if such an omnivory link 
stabilizes the system by reducing the perturbations induced by shifts in death rates. The shortcut 
link modifies the rate equations for the u level by adding a gain term 

OL")(OjX w X u 

and for the w level by adding a loss term 

As in the case of self-interaction, the dimensionless factor a defines the intensity of a shortcut 
relative to regular predation. The gain term has a conversion factor x which is distinct from A 
since preying across several trophic level can have different efficiency in terms of the birth rate of 
the predator, compared to preying on the next nearest level below. The uth and wth equations for 
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the stationary concentrations are modified accordingly, 

d _ 

a 

Ax„,_i = x w+1 H h ax u (41) 

a 

As in the case of self-interactions, it follows that the omnivory link between u and w directly 
affects the population of the levels that are recursively downstream from u and w. Similarly to 
the self-interaction case, we denote corrections to the population steady states resulting from the 
effect of a death rate shift on a food chain with shortcuts by Ax fc . A formal recursion equation for 
Ax k , analogous to Eq. (|30l ) reads 

AAx fe _i = Ax k+1 + 4, P e + S k>w aAx u - 8 k>u xaAx w (42) 

Since the perturbed level, the upper level of a shortcut, and the lower level of the shortcut can 
belong to either {n} or {n — 1} partition, there exist 2 3 = 8 distinct scenarios. Below we present 
detailed analysis for two examples of enhancement and dampening of perturbations by a shortcut. 
The analysis of the remaining 6 scenarios is straightforward and yields additional stabilizing and 
destabilizing examples. 



A. Example of destabilization 

We consider a case where all 3 levels p, u, w belong to the {n} partition, as illustrated in Fig.[U 
Let us first assume that the death rate shift occurred at p > u. The corrections to population 
induced by the omnivory link depend on Ax u and Ax w , which, in their turn, could be expressed 
through x w -i via the relation (|26|) linking both partitions. 

Ax w = -^V\Ax w -i, Ax u = \ {u ~ w)/2 Ax w . (43) 

On the other hand, Ax m _i depends on Ax w directly (see Eq. (|42l) for k = w), and on Ax u _i 
recurrently through Ax w+X . Ax n _i in turn directly depends on Ax u via Eq. (l42l) for k = u, and 
recurrently through Dx u+ i, which depends on the death rate shift at p. Putting together all these 
dependencies yields an equation for Ax w _i and Ax w , 

Ax w ^ = eX^- 2 ^ 2 + (A(«^)/ 2 - xA^")/ 2 ) . (44) 
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FIG. 8: The effect of omnivory interaction between levels u = 3 and w = 1 on the propagation of a 
perturbation caused by a shift by ea of the death rate at level p = 5; all these 3 levels {u, w,p} are in {n} 
partition. 

Combining Eqs. (I43ll44l) produces an expression for the population correction i < w, 

\(i-p-l)/2 

At. = (45) 

1 1 + a* (A(«-™)/ 2 - x \ {w - u)/2 ) ' 



eA^'-f- 1 )/ 2 

AXj = ~ 7 l + a * (A(«— )/ 2 - X A(—")/ 2 ) ' (46) 



and, via (1261) , for any j, 



For w < i < u, only the u and p levels affect the population perturbation, so that, 



Axi = ttt, (47) 

1 - a*x\ [w ~ u)/2 

and for i > u no effect of the omnivory link is present. These results hold for arbitrary positions 
of the death-rate-shifted level p with respect to u and w: Once p G {n}, all {n} partition levels, 
including u and w, are perturbed. The only difference with the above case occurs when p < i. 
In this case, the direct effect of death-rate perturbation eA^ _p_1 ^ 2 should be subtracted from Ax^ 
defined via Eqs. (I45ll47l) . It follows from Eqs. (I45ll46l) that the sign of the effect depends on the 
values of A and \- I n a common scenario with A < 1 and \ ~ \ me effect of the omnivory 
link at the i < w and all j levels is destabilizing: the perturbation induced by a shift in death 
rate is enhanced in the presence of the omnivory link. However the strongest enhancement or 
destabilization, which is universal and manifests itself for arbitrary A and x> occurs for the w < 
i < u levels of the {n — 1} partition, Eq. (PTTl) . In the idealized case with A = % = 1, the two terms 
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FIG. 9: The effect of the omnivory interaction between levels u = 3 of the {n} partition and w = of the 
{n — 1} partition on propagation of perturbation caused by a shift by ea of the death rate at level p = 5, 
which belongs to the {n} partition. 

in brackets in the denominators of (1451 l46l) cancel each other, so that the destabilizing effect of the 
omnivory can be seen only at the levels with w < i < u. 



B. Example of stabilization 

Here we consider a case where the omnivory link connects a level u of the {n} partition to a 
level w of the {n — 1} partition, while the death rate perturbation occurs at a level p of the {n} 
partition, as illustrated in Fig. |9j First consider the case when p > u. From (1421) it follows that 

Ax n _! = eX aj - p - 2)/2 - Ax w a X /X (48) 

Taking into account that Ax w is related to Ax u ^i via a simple recurrence, 

Ax w = Ax n _iA (w ~ n+1)/2 , 

one obtains for i < u 

AXl = 1 + a X ^ w - u -^ 2 (49) 



Consequently, for j < w the perturbation is a simple recurrent continuation of (1491) . 

A Xj = ^— (50) 

For Ax^+i it follows from (|4T|) that 



AAx^_! = Ax w+1 + aAa; u = Ax w+1 (l + a\ {u ~ w ~ l)/2 ). (51) 
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Hence for j > w 

At ■ = ' (52) 

3 (1 + a X X^- u -^ 2 ) (1 + aA(«— !)/ 2 ) ' ' 

A similar scenario occurs when w < p < u; to get the population corrections one needs to subtract 
the direct effect of the death rate shift eA( l-p ~ 1 " 2 from the corrections Axt given by Eq. (|49l) for 
i > p. When p < w, there is no effect of the shortcut on the propagation of perturbation in 
the {n — 1} partition and the lower part, j < w, of the {n} partition. However, for j > w the 
shortcut does have an effect on perturbation amplitude, which is expressed by the second term in 
the denominator of (l52l) . 

e7 ^_xi-p-i)/2 

Ax-j = —- r? ttat (53) 

Thus, we observe that an omnivory link between u E {n} and w E {n — 1} stabilizes the 
food chain by reducing the effect of a death rate shift on the population. This conclusion holds for 
arbitrary A and x> and for any relative position of the levels linked by omnivory and the level at 
which the death rate perturbation occurs. 

VII. DISCUSSION AND CONCLUSIONS 

Using linear Lotka-Volterra food chains, we have investigated how perturbation at one trophic 
level cascades to other trophic levels, and how such perturbation cascades are affected by the 
introduction of additional links in the food chain, representing either self-interaction (intra- specific 
competition or cannibalism), or omnivory. Our results show that such additional links can either 
dampen or enhance the perturbation cascade, depending on the position of the level at which the 
perturbation occurs, and the levels connected by the additional link. Thus, self-interaction and 
omnivory can be either stabilizing or destabilizing, depending on the details of the setup. 

In establishing these results, we obtained a closed-form solution for the steady state populations 
of Lotka-Volterra food chains of arbitrary length, and we derived analytic expressions for the 
changes in the steady states in response to shifts of death rate, self-interaction and additional 
intraguild trophic relations. 

Our definition of stability, based on susceptibility of steady state populations to long-term re- 
moval or addition of a population of a certain level is somewhat different from previously con- 
sidered ones. These previously considered definitions of stability include the size of cascading 
extinction events following the removal of randomly chosen species [6|], permanency, or local- 
ization to a certain region of phase space for systems with possibly non- stationary asymptotic 
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behaviour Hi 111 . time of relaxation to a new equilibrium [50, and the mere convergence to a steady 
state 1120. Our definition is inspired by practical applications in which the effect of continuous 
harvesting of a certain species on the whole food web needs to be considered. However, suscep- 
tibility to perturbation is intrinsically related to other definitions of stability, such as the existence 
of steady states and lack of extinction events. In response to a perturbation, a less susceptible (and 
thus more stable, according to our definition) system deviates less from a steady state than a more 
susceptible one, thus having a lower probability of leaving the basin of attraction of a stable fixed 
point, or cross the boundary of a phase region with permanent behaviour. This can be illustrated 
by considering, for example the expressions for stabilizing (1331) and destabilizing (|37l) effects of 
self-interaction. When Xj + Axj is negative, a steady state is no longer possible, which means 
that either the population of the level j(and consequently, all higher levels) becomes extinct, or the 
behaviour of the system becomes non- stationary. When the effect of additional links is stabiliz- 
ing, (|33l) . this happens for a higher level of harvesting, and when the effect is destabilizing, (|37|) . 
the sustainable level of harvesting is lower. In other words, if a system is more stable according 
to our susceptibility criteria, it would also be more stable according to the traditional definitions. 
Hence our conclusions about stabilizing (or dampening) and destabilizing effects of competition 
and intraguild predations are relevant for other definitions of stability. 

For our analysis, we considered topologically very simple systems of a linear food chains with a 
single additional link representing intraspecific competition or omnivory. In addition, we assumed 
linear functional responses to describe the predator-prey relation between adjacent trophic levels 
in the food chain. Nevertheless, our results can be generalized in a number of ways. 

For example, numerical simulations reveal that our results remain true qualitatively with non- 
linear functional responses, as is shown in the Appendix. In Fig. [TT] we present an example of a 
destabilizing effect of an omnivory link, and in Fig.[[2]an example of stabilizing effect of omnivory 
is shown. It follows from these figures that in systems with Type-II functional response the effect 
of omnivory links on stability is qualitatively similar to the structurally equivalent Lotka-Volterra 
counterparts. Moreover, our results can be generalized to more complicated forms of regular 
food webs, such as the Cayley trees illustrated in Fig. [lOl In "upwards" Cayley trees (Fig.fTOla). 
m predators in the next higher trophic level are feeding on the same prey, and in "downwards" 
Cayley trees (Fig.[K)]b), a given predator feeds on m prey species on the next lower trophic level, 
so that for m = 1 the linear food chains are recovered. 

For example, for upwards Cayley trees (Fig. [10] a), the population dynamics at the various 
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FIG. 10: Sketch of branching food webs (Cayley trees) in which \i = 2 predators are feeding on the same 
prey, a), and in which each predator is feeding on ji = 3 preys, b). 



trophic levels is given by 



dxk 
~dt 



x k (\ax k -i - fJ,ax k+ i - 8 k ) • 



(54) 



In this case, the expression CCD) for the steady state populations still holds, but in all subsequent 
steps of the analysis presented in Section II, A should be replaced by X//i and d by d/fj,. In the 
complimentary scenario of "downwards" branching (Fig.fTOb). A in (flTT) and subsequent equations 
should be replaced by //A. Starting with these expressions, it is possible to derive results for the 
effects of self-interaction and omnivory on food web stability in Cayley trees, that are analogous to 
the results presented in Section III - VI. In particular, omnivory and self-interaction can have both 
stabilizing and destabilizing effects in these more complicated food webs. Extrapolating from this, 
we conjecture that in more realistic food webs with more complex topologies and with a multitude 
of self-interaction and omnivory links, it would be even harder to draw general conclusions about 
whether a particular link stabilizes or destabilizes the whole food web. 

Our results may help to resolve the long-standing controversy on the effect of omnivory and 
self-interaction on the stability of food webs by essentially showing that general uniform con- 



clusions regarding stabilization l^l jji fioi 



121] or destabilzation 



14|1 cannot be drawn. Rather, 



whether a particular link in a food web leads to more or less stability very much depends on the 
details of how this link is embedded in the whole food web. For example, the sign of the effect 
depends on such features as the relative position of omnivory or self-interaction and the perturba- 
tion level from the top of the food web. In our simplified modeling situation, these positions are 
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unambiguously defined. However, in reality, the top predator is often rare (e.g. big carnivores, 
which may often have low population densities), which will make the determination of positions 
relative to the top level difficult in real food webs. Nevertheless, it seems feasible that our results 
could be tested experimentally in simple and well-compartmentalized trophic systems. 



VIII. APPENDIX 



A. Type-II functional response 



Here we show that the results derived above for the Lotka-Volterra type of functional response 
qualitatively hold for a more realistic, Type-II functional response. We consider a linear food 
chain with nearest neighbour predator-prey interactions so that a rate of change of population 
density Xk(t) of the kth trophic level is given by 



dt 



■J'k 



0&+i,fc£fc+i 



— 5i 



(55) 



Here a it j are predation strength of species i on species j, is the conversion efficiency which 
connects the birth rate for the species k to the amount of other species it consumes, h). is the 
handling time for species k, and d& is the death rate coefficient of the species k. Evidently, the 
top predator, occupying the nth level, does not have any species that prey on it, which could be 
expressed by setting x n+1 = In addition, the basal species, occupying the basal trophic level, 
is characterized by the logistic growth term, which mimics the finite input of the energy into the 
system (or finite carrying capacity K). The linear death term d is absorbed into the linear part of 
the birth term f3, which gives for the rate of change of the basal species concentration, 



dx 
~~dt 



x 



a\nXi 



(56) 



1 + a lfi h 1 x _ 

As in the main text, we consider the effect of "shortcuts", or omnivory food links, reflecting 
predation of the species u not only on its nearest neighbour u — 1 below it in the food chain, but also 
on the species w with the intensity a and conversion efficiency \ (which replaces A). Introduction 
of such link modifies the rate equations (155156b for the species w, 



dt 



dec UJ 
~dt 



i + h 



(57) 
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FIG. 1 1 : The ratio lu^ 



K- 



•a;*. J 



of the differences between the perturbed and unperturbed population 



with (a = 0.2) and without (a = 0) omnivory, for Lotka-Volterra functional response (circles) and Type-II 
functional response with h = 0.25 and x = A = 1 (squares) and \ = 0-5, A = 1 (diamonds). 



the species u — 1 



u-1 



u-1 



eft 



"I - Q"u,u— l%u%u—l 



1 



1 "I - Q-m,u— \h u x u — i 
1 



and for the species w 



~dt 



1 + h u (aax w + a niU _ix u _i] 

(^U,U— 1 ^M^ll- 1 



(58) 



1 H - ^u.u—l^u^u- 



+ 



•\i^U,U— l-^u— 1 "I - X^^-^io 
1 -(- ^iufciu.u— l*^u— 1 "I - Cy(ZX„ 



(59) 



We use the following values of the constants, a 



ij 



1, \ = 0.75, 1, Xi = 0.5, 0.75, 1, 



K = 1, P = 1, 5i = 0.1, and hi = 0.25 for all i and j. As an example of destabililizing effect of 
intraguild predation, Fig.QTJ we consider a linear food chain with an omnivory link of the intensity 
a = .2 between levels u = 3 and w = 1. The perturbation is caused by a shift by ea, e = 0.1 of 
the death rate of the level p = 5; all these 3 levels {u, w, p] are in {n} partition. 

As an example of stabilizing effect of intraguild predation, Fig. [121 we consider an omnivory 
link between levels u = 3 of the {n} partition and w = of the {n — 1} partition. The perturbation 
is caused by a shift by ea, e = 0.1 of the death rate of the level p = 5; which is in {n} partition. 

It follows from Figs. [TH [12] that the conclusions made in Section [VI] qualitatively hold for 
foodwebs with Type-II functional response as well: Depending on the position of the level with a 
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FIG. 12: The ratio uj^ = ^hhwm Q f the differences between the perturbed and unperturbed population 
with (a = 0.1) and without (a = 0) omnivory, for Lotka-Volterra functional response (circles) and Type-II 
functional response with h = 0.25 and x = A = 1 (squares) and \ — I = 0.25 (diamonds). 

death rate shift and a shortcut link, the omnivory can either stabilize of destabilize the food web. 
When the omnivory link connects the levels u = 3 and w — 1 of the {n} partition and the death 
rate shift occurs at the level p = 5 also belonging to the {n} partition, the perturbation is enhanced, 
especially at the level(s) w < i < u of the {n — 1} partition. Conversely, when the omnivory link 
connects the levels u = 3 of the {n} partition and w = which belongs to {n — 1} partition 
and the shift of the death rate occurs at the level p = 5 which belongs to the {n} partition, the 
perturbation is dampened. 
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